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ABSTRACT 

The association of all three soft gamma-ray repeaters (SGRs) with supernova remnants 
has established that SGRs are young neutron stars, and has given us a starting point 
for detailed modeling. One of the most popular classes of models involves strongly 
magnetised neutron stars, with surface dipole fields B ~ 10 14 — 10 15 G. In such 
strong magnetic fields, many otherwise negligible processes can play an important role. 
Here we consider the effects of vacuum polarisation on Compton scattering. Vacuum 
polarisation introduces a characteristic density-dependent photon frequency at which 
the normal modes of polarisation become nonorthogonal and the mean free path of 
photons decreases sharply. Our analytic results and Monte Carlo simulations of photon 
propagation through a magnetised plasma show that this effect leads, under a wide 
range of physical conditions, to a broad absorption-like feature in the energy range 
~5 keV— 40 keV. We discuss this effect in light of the spectra from SGR 1806-20. 
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1 INTRODUCTION 



The first soft g amma-ray repeater (SGR) to be discovered, 



SGR 0525-66 (Mazets fc Golenetskii 1981; 



Mazets et al 



|l97S| ), was shown soon after its discovery to be positionally 
coincident with th e N49 supernova re mnant in the Large 
Magellanic Cloud (Evans et al. 1980). More recently, the 



despite the highly super-Eddington luminosities, the plasma 
is probably confined near the surface of the neutron star. 
Magnetic fields could confine the plasma, and perhaps sup- 
press plasma instabilities, with pulsar-like field strengths of 
~ 10 12 G (Lamb 1982; Katz 1982, 1993, 1994). 

There are lines of argument suggesting much greater 



other two SGRs (SGR 1806-20 and SGR 1 900+14) have also 



been associated with supernov a remnants (Kouveliotou et al 
1991; iKouveliotou et al. 199 [Hurley et al. 1994 |Kulkarni 



field strengths, on the order of 10 14 
in Thompson and Duncan (1995), a <■ 



- 10 b G. As suggested 
6 x 10 14 G field could, 



&i Frail 1993; |Kulkarni et al. 1994[ ). These associations have 
led to a consensus that SGRs originate from young neutron 
stars. 



The peak luminosities inferred for bursts from SGRs are 
typically 10 40 — 10 42 erg s _1 for ~ 0.1 — 1 seconds (except- 
ing the unusual March 5, 1979 event from SGR 0525-66, 
which had an inferred peak luminosity of ~ 2 x 10 45 erg 
s~ ). These luminosities are ~ 10 2 - 10 4 times the Edding- 
ton luminosity L E = 1.3 x 1O 38 (M/M ) erg s" 1 at which 
the radial radiation force from Thomson scattering exactly 
balances the gravitational force on a fully ionised hydro- 
gen plasma. Nonetheless, the spectrum from SGR 1806-20, 
which has burst more than 110 times and is thus the best- 
studied of the SGRs, appears not to undergo the reddening 
one wo uld expect if the burst created an expanding opti- 



cally tlick fireball ( Fenimore, Laros, fc Ulmer 1994[). Thus, 



within the ~ 10 year lifetime of the N49 supernova rem- 
nant, spin a neutron star down to the 8 s period seen in 
the March 5 event. Moreover, a superstrong magnetic field 
may be able to account for the energetics and time scales of 
SGRs (see Thompson & Duncan 1995 for one model). These 
models are extremely complicated, and there are many as- 
pects of microscopic physics that are important only in very 
strong magnetic fields. Here we concentrate on the effects 
of vacuum polarisation and, in particular, on the spectral 
effects of the second vacuum frequency. 

2 3 

In strong magnetic fields B>B C , where B c = —i— = 
4.414 x 10 13 G is the magnetic field at which the electron 
cyclotron energy equals the electron rest mass energy, vac- 
uum polarisation from virtual e + e~ pairs can become sig- 
nificant (see, e.g., Meszaros 1992 for a discussion). In fact, 
for strong enough fields the vacuum contribution to the di- 
electric tensor can exceed the plasma contribution ( Pavlov 
fc Shibanov 1979|; [Ventura, Nagel, fc Meszaros 1979^ Pavlov 
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& Gnedin 1984). This has effects on the normal modes and 
cross sections for photons. Well below the electron cyclotron 
frequency ujb = eB/m e c, the polarisation normal modes for 
photons in a very strong magnetic field are close to linear 
over most photon frequencies and incident angles. 

At most photon frequencies uj <C ujb, the scattering 
cross section depends strongly on the polarisation of the pho- 
ton. If the electric field vector of the photon is in the plane 
formed by the magnetic field and the photon propagation 
direction k, the photon is in the parallel mode and the scat- 
tering cross section is <ti ~ or, where <tt = 6.65 x 10~ 25 cm 2 
is the Thomson cross section. If the electric field vector of 
the photon is perpendicular to the B — k plane, the photon 
is in the perpendicular mode and 02 ~ {lu/lub)' 2 ot -C ot', 
this reduction in cross section is understandable, because in 
a strong magnetic field it is easier to oscillate an electron 
along the field than to oscillate it across the field. 

There is, however, a small range of photon frequencies 
in which the vacuum and plasma contributions to the di- 
electric tensor are comparable to each other and the normal 
modes become strongly nonorthogonal over a broad range of 
angles. This frequency, called the second vacuum frequency 
(the first vacuum frequency is near the cyclotron frequency) , 
depends on the magnetic field B and plasma number den- 
sity N: uj c2 ~ iV 1 / 2 ^- 1 for B < B c and uj c2 ~ N 1/2 B' 1/2 
for B 3> B c . Vacuum polarisation does not change the 
polarisation-averaged cross section, but it does change the 
normal modes and redistribute the total cross section be- 
tween them: when u> ~ uj c2 , o\ ~ a 2 ~ o~t- At any particu- 
lar density the cross section is unusually high for a narrow 
range of frequencies (see Figure 1). Conversely, if we assume 
that most of the flux is transported in the lower cross sec- 
tion mode, for any photon frequency there is a density at 
which the cross section is high. Thus, the optical depth seen 
by a photon is greater than it would be in the absence of 
the enhanced cross section at uj c2 . Because the resonance 
frequency varies with density, in an atmosphere of varying 
density the vacuum resonance has an effect on the spectrum 
in a wide range of photon energies. 

Here we expand on results first presented in Miller & 
Bulik (1996). In § 2 we give the normal modes and cross sec- 
tions for photons in a strongly magnetised plasma, includ- 
ing the effects of a significant electron-positron pair density 
and of magnetic fields B > B c . In § 3 we give a qualitative 
derivation of the upper and lower bounds to the vacuum res- 
onance feature as a function of magnetic field, scale height, 
and the fraction of positrons. In § 4 we describe our numer- 
ical method for following photons as they scatter through 
the atmosphere. We also give the results of these simula- 
tions, and show that an absorption-like feature is present at 
reasonable optical depths for any magnetic field B>B C . We 
neglect the effects of photon splitting; at the energies Tiu)< 
50 keV of interest here the probability of photon splitting 
within ~ 10 6 cm is small (see Adler 1971). We likewise ig- 
nore proton cyclotron scattering, because we focus on scale 
heights Z< 1000 cm at a constant magnetic latitude, imply- 
ing that proton cyclotron scattering will yield extremely nar- 
row features, , which would be unresolvable by current 
detectors. Finally, in § 5 we discuss how our results, com- 
bined with observations, can constrain physical parameters 
in SGRs. 



2 VACUUM POLARISATION EFFECTS 

It is convenient to describe radiative transfer in a magne- 
tised medium in terms of normal modes. Including both the 
electric permittivity tensor e and the magnetic permeability 
tensor fj,, the dispersion equation becomes 

UJ 



k x [/T 1 (k x E)] + -eE = 0, 



(1) 



where E is the electric field. A simple and elegant way 
of solving equation (|l|), describing the normal modes, 
and calculating cross sections has been introduced by 
Gnedin & Pavlov (1974). Given a system of coordinates with 
the magnetic field along the z-axis, and the wave vector in 
the yz plane at an angle 9 to the magnetic field, we define 



1 + i (dcx + e vv cos 2 9 + e zz sin 2 1 



— 1 _ —1 2 

~f^xx L^yy COS ( 



III 



and 



nc = - (e xy cos 9 + e xz sin ( 



fi zz sin 2 ( 



: sin 26> 



fi xz sin 26») 



, sin 29 



fizz sin 2 9 + u xz sin 29) 



(2) 



(3) 



(4) 



where fi a £ is the ab component of the inverse fi tensor. The 
roots of the dispersion equation are 



= m± yJn 2 L +n\ 



(5) 



where j — 1,2 indicates which normal mode is selected. 
The refraction coefficients are given by the real part of rij, 
and the total absorption coefficients are proportional to the 
imaginary part of nj \ £j — (2u)/c)lm(nj), and the cross sec- 
tions in the two modes are 



<jj — (2uj / cN e )lm(iij) . 



(6) 



The polarisation of the normal modes can be described by 
the position angle \j between the major axis of the polari- 
sation ellipse and the projection of the magnetic field B on 
the plane perpendicular to the wave vector k, and the ellip- 
ticity V , whose modulus is equal to the ratio of the minor 
axis to the major axis of the polarisation ellipse, and whose 
sign determines the direction of rotation of the electric field. 
The ellipticity and position angle are 



r 7 + 1 



r,- exp(2ixj) 



-nc 



±y/n 2 L +n: 



til 



(7) 



The polarisation vectors can be described by a complex vari- 
able 6, or by two real parameters q and p 

b = q + ip^^±. (8) 

nc 

A detailed discussion of the polarisation of the normal modes 
is given by Pavlov, Shibanov, & Yakovlev (1980). 

Using the cyclic coordinates e± = 2~ 1 / 2 (e x ±ie y ), eo = 
e z (Meszaros 1992, p. 95), the polarisation vectors are 



ei = ( ' . I\ >in (I 
where 

Kj = b\i + (-iy (1 



hCje*** (K jCOS e±l) 



-2\l/2* 



(9) 



(10) 
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and Cj = (1 + j-Kj| 2 ) 1//2 is the normalisation. It should be 
noted that in general Kj is complex. 

2.1 Vacuum effects in the absence of matter 

In the presence of a strong magnetic field the vacuum ac- 
quires a non-zero polarisability tensor due to the existence 
of virtual electron positron pairs, so that the permittivity 
tensor etj and the magnetic permeability tensor fiij of the 
vacuum become different from unity. Charge symmetry im- 
plies that both tij and fXij are diagonal in a vacuum. 

For example, if B <C B c and huj <C mc 2 , so that photon 
splitting and the creation of real pairs can be neglected, the 
dielectric tensor and magnetic permeability tensor are given 
by 



and Li = | + /„* dx lnTi(l + x) « 0.249. The asymptotic 
limits of the functions r)\ t 2{h) have been found by Tsai & 
Erber (1975). In the low field limit, B « B c , 



14, 2 13 4 

_8_, 2 _ 379 l4 
45 5040 ' 



and when B ^> B, 
2 



Vi(h) 



-ft- 



(21) 
(22) 

(23) 
(24) 



m{h) » --ft^ln^ft). 

Using equations ( |l4| ) and ( |l5|) we can rewrite equa- 
tions (|l3|) 



-1 -, 
Ma6 - 1 



•2a 










-2a 


») 








5a / 


•2a 











-2a 











—6a 



(11) 



(12) 



a ( B \ 2 , . , „ 

where a = — — and a is the fane structure constant. 

457T \BqJ „ „ „ 



Inserting equation (\12\) into equations 
obtain 



and 



71! 



vac 
vac 

n c 



1 4 a sin 9 . 

4 



3 . 2 
—asm 9 , 



(13) 



0. 



vac - 1 - / iiac , vac\ 

nj = g C n i + n 2 ) 

= l+^sin 2 0(7 ? i(/i) + 772(h)) (25) 

nl ac = i(nr c -n2 ac ) = ^- sin 2 6»(r?i (ft) - 772 (ft)) (26) 

<? c = . (27) 

The normal modes are polarised linearly for all angles of 
propagation except for the degenerate case of propagation 
along the direction of the field when they are circular. Such a 
behaviour could also be found by considering the symmetry 
of the problem. 



2.2 Propagation of radiation in a system with 
matter present 



The roots of the dispersion equation in a vacuum are real, 
and the refraction indices of the normal modes are 



vac vac . vac i , r» 2 , 

Tlx = n i + n L = 1 + 2a sm l 



vac vac 

n 2 = nj 



-u a c 



= 1 + ^asin 2 ( 



(14) 
(15) 



Note that knowledge of the refraction indices n" ac and rC^ ac 
suffices to determine the influence of vacuum polarisation 
on the propagation of light in the presence of matter. 

The refraction indices nl°2 have been calculated for 
hui <C m e c 2 and arbitrary magnetic field by Tsai & Erber 
(1975), 



vac i . / 7 \ . 2 n 

rii = 1+ -r- 771 (ft) sin 9, 

47T 

vac t . & / 1\ • 2 /1 

n 2 = 1 + — 772 (ft) sin 9, 

47T 



(16) 
(17) 



where ft • 

mW = 



B_ 

K 



, and 



81nri(l + ft)-4ftlnr(l + ft)- -*(l + ft) 



Consider now a system in which the plasma density is low 
enough that the contribution of plasma to the dielectric ten- 
sor is |eL — 1| <C 1, and the magnetic field is small so that 
|e|/ c - 1| < 1 and I//7. 1 - 1| < 1. This condition is satisfied 

B , 
for — < 40 and u p <C ujblj, which holds for 

Be 



(t)- 



(28) 



^»l-6xlO- 8 ( T - i ^ 
mc 2 WO^cm- 6 

Such a system can be described by the dielectric tensor 

tab = <="aT + tab ~ &ab , (29) 

and the vacuum permeability tensor. The plasma dielectric 
tensor can be expressed as e a t = 5 a b — [ — S ] fl a 5, where 



{47vNe 2 Y /2 

u)„ = is the plasma frequency and U a b is the 

V m / 

plasma polarisation tensor. The plasma polarisation tensor 
is diagonal in the cyclic coordinates, and for cold electron 
plasma is given by 



-2ft In ft - 2ft 2 + 2ft ln(27r) + — + - - 8L1 (18) 



??2(ft) = -4ftlnr(l + ft)+4ft 2 >I'(l + ft)+2ftlnft 
+2ft[ln(27r) - 1] - 4ft 2 + | . 



In these expressions $(a; 

fx 

InF^x) = 



dt lnr(t) + -(x - 1 - ln(27r)), 



(19) 



(20) 



I I n r 



(30) 



lu 4- olujb — Z7r 4- au>B 

where a — — 1,0, +1, and j r = (2/3)(e 2 /mc 3 )ti; 2 is the ra- 
diative width, and we denote to t = uj — i'yr- When there 
is an admixture of pairs, such that the positron fraction is 
/ = n c + /(n e - + n c +), the cyclic components of the polari- 
sation tensor are modified 



nL = /n Qa + (i - /)n_ 



(31) 
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Inserting equation 
obtain 



into equations 



and 



and 



nc 



1 + 



1 u v 



(1 + cos 2 9)- 



4 to 

+ ^sm 2 6(r ll (h)+ri2(h)), 



i • 2 a 1 

+ sin 6* — 



"2 w 



^72 W) 



(32) 
(33) 



(34) 



When / 7^ 0, the parameter is not changed, whereas 



V 2 iw 2 , 



U)U)B 



■ COS 6> . 



nc 

Inserting equations (^) and (|3^) into equation 
giecting terms proportional to 7 2 we obtain 

1 „:_2, 



(35) 
and ne- 



q = 



and 



1-2/ 2 cos 6» w 



X 



2tv Up 



P = 



LUB'fr 



1 - 2/ 2 cos 6> lj 2 



(36) 



(37) 



The normal modes are linear and orthogonal for most fre- 
quencies and angles, since q 2 ^> p 2 for most of the parameter 
space. The absorption coefficients in this regime can be ob- 
tained by inserting equations (^2|), (^3|),and ( |34| ) into (JEJ), 
and noting that the ratio ^ is large and real. Thus, we 
obtain: 



£2 



= Nar sin 6 
= N c t ( — 

\0Jb 



(38) 
(39) 



The critical points are where q(to) = 0. Such points exist 



when 

Q w s r cm 
5^7? fa CO 



>1. 



(40) 



Generally, for a system consisting of electron plasma and 
magnetised vacuum there exist two frequencies for which 
q(u>) — 0, one very close to the cyclotron frequency, and the 
other at u> <C cub- In this work we will be interested in the 
latter. At this point, called the second vacuum frequency, the 
normal modes are circular for most propagation angles. This 
leads to an increase of the absorption of the low cross section 
mode and a simultaneous decrease of the absorption in the 
high cross section mode. From equation (pS6|), the second 
vacuum frequency is 



0J C 2 = Wp 



V2(h)) 



-1/2 



with the limits 

hco c2 = fiuj v (15tt/ a) 1 /2 (B c /B) 
= 30N^ 2 (Bc/B) keV 



(41) 



(42) 



for B < B c and 

ftw c2 = H (37r/a) 1/2 (B c /B) 1/2 



13iV2 6 /2 (B c /B) i/2 keV 



1/2 



(43) 



for B > B e , where iV 2 6 = A^/IO 26 cm -3 . Polarisation of the 
normal modes at uj C 2 is determined by the parameter p. The 
modes are circular when p 2 < 1 and linear when p 2 > 1. 
In the case when p = 1 and q = the two polarisation 
modes coincide. Inserting equation ( ffl| ) into (|37 ) we obtain 
the value of p at the vacuum frequency: 

™s 1 



jj(u c2 , 0) = 



(44) 



1 - 2/ 3 cos mc 2 ' 

We can define a critical angle 9 C as an angle for which 
p(u! C 2,9 c ) — 1. For the magnetic fields considered here 8 C 
is very close to n/2; n/2 — 8 C m (a/3VB/B c ) for elec- 
tron plasma. Again inserting equations (pS2|), (pS3|), and ( p4"| ) 
into (g), and noting that the ratio ^ is imaginary and small 
for the angles 9 < 9 C , we obtain 



iV(TT 




.(45) 



In the case of propagation orthogonal to the field, when 
9 > 9 C , the normal modes become linear and the absorption 
coefficients are given by equations ( ^ ) and d3s|). 

The width of the resonance is associated with the extent 
of the interval where q(u>, 9) < 1. It can be estimated as the 
inverse of the derivative of q at the vacuum frequency 

-1 * 2 

w c2 



dq 
dui 



(1-2/)- 



LOB 



(46) 



The vacuum frequency and the shape of the vacuum reso- 
nance do not depend on temperature since thermal effects 
alter the dielectric tensor weakly for uj <C ujb- Unlike the cy- 
clotron resonance, the vacuum resonance is a result of global 
properties of the medium and does not depend on any indi- 
vidual particle resonance. 

2.3 Scattering cross section 

Let us now consider scattering in a hot, magnetised plasma. 
A photon with energy u), wave vector k propagating at an 
angle 9 to the magnetic field B, and polarisation vector e 
scatters on an electron with parallel momentum p. After 
the scattering the photon has energy u>' , wave vector k' , 
polarisation vector e', and the electron momentum is p' . 

The scattering cross section is given by (Meszaros 1992, 
p. 131) 



da(u>, 9,p>—> u>' , 9' ,p 



2 U> 



u) \k- ft' I 



|(e|n|e')| 



(47) 



where ro = e 2 /mc 2 is the classical electron radius. After 
averaging equation (^) over <f>' we obtain 



da (ui , 9 , p h- > u> ,9 ,p ) 



2 lu m 



X 



ui \k — k'\ 

x^|n ±tt | 2 |e Q | 2 |eL| 2 , (48) 
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Figure 1. Left panel: Scattering cross section, in units of the Thomson cross section, versus photon energy. The solid lines correspond 
to the total cross sections in the two modes. The dotted line represents the cross section for the low mode to the low mode, the long 
dashed line represents the low mode to high mode cross section, and the short dashed line represents the high mode to low mode cross 
section; the high to high cross section is essentially equal to the total high mode cross section. Note that the width of the resonance in 
the low to high mode cross section exceeds that in the low to low cross section. The inset shows a close up of the frequency range near 
the resonance. Right panel: Mode switching probabilities for the low mode (solid line), and the high mode (dashed line). 



where we have the + sign for scattering on electrons and — 
for scattering on positrons. The scattering amplitudes, with 
first order relativistic corrections, are 



n+ = l- 
n- = f + 



LOB 



u' +UB- pk' + k' 2 /2 + i-y r 
l£B 

U) — lob — pk — fc 2 /2 + ry r ' 



(49) 



f + 



(p + k/2)(p-k' + fc/2) 

LO — pk — fc 2 /2 + 47r 

(p-k'/2)(p-k' + k/2) 
to 1 - pk' + k' 2 /2 + i-) r 



Assuming coherent scattering, i.e. k = k! , we average over 
the final states and use the polarisation normal modes to 
obtain the total cross section from one mode to another: 



© 0000 RAS, MNRAS 000, 000-000 



6 Bulik and Miller 



|n±o; | . 



(50) 



where A J a = | J d cos 9\e 3 a \ 2 . In general the integrals A 1 a 
are difficult to calculate analytically. We can introduce the 
mode switching probability for the two modes 



Pn = 



+ a j 



(51) 



We will find the cross section and the mode switching prob- 
abilities in some limiting cases. 

When polarisation modes are linear for most of the an- 
gles, e.g., outside the resonance, where polarisation is vac- 
uum dominated, we have \b\ — \riL/n c \ ^> 1, we find using 
equation @ that Ki -> 26, K 2 -> - I /2b, C x -> I /2b, and 
C?2 — ► 1. The components of the two polarisation modes are 
then 



l 



and 



2 



V2 



cosSe^"^ , el ~ sin( 



V2 



(52) 



(53) 



The total scattering cross sections in the limit u = 
(ju) S 5, 1 °f photons at initial incident angle 9 to the 
magnetic field are 



&1-.2 



C2—1 



(Tt sin f , 

3 1 2, 
— (Tt— cos l 

4 u 

1 1 

-<TT- , 

4 u 

3 1 
-or— • 

4 u 



(54) 
(55) 

(56) 

(57) 



Thus the probability of switching from the high cross section 
mode to the low cross section mode is 



■ cos 2 9 



Pl-,2 



u sin 2 9 + | cos 2 I 



3 / u 
4Wb 



cot 



(58) 



where the last approximation holds when 9 ^> \f\-/j^- The 
probability of switching from the low cross section mode to 
the high cross section mode is 

1 



P2-.1 = 



(59) 



In contrast, when 6 « 1, which holds for most of the 
propagation angles in the resonance (9 < 9 C ), the modes 
are circular, and Ky — > 1, K 2 — > — 1, C\ — > l/y/2, and 
C?2 — ► 1/ \/2. The components of the two modes are then 



i 



and 



2 



i^e^'^l =F cosfl) , » ^= sinf 

2 V2 



±-e 



1 ± cos 6*) , e z 



^2 



sin 9 



(60) 



(61) 



The total scattering cross section in the limit u 3> 1 for 
any initial or final polarisation state is then, in the limit 



sin 9 ^> u 



-1/2 



1 • 2, 

— ctt sin l 

4 



<Ti->f 

The mode switching probabilities are 



(62) 



Pi^2 = P2^i = ^ • ( 63 ) 

We calculate the polarisation of the normal modes using 
the cold plasma approximation, since the thermal effects are 
not important in the range of energies far from the cyclotron 
resonance. We then use equation ( ft8j ) and the scattering am- 
plitudes ( |49| ) to find the cross section as a function of the 
final angle. We present the scattering cross section as a func- 
tion of photon energy in the left panel of Figure 1. In the 
right panel of Figure 1 we show the mode switching prob- 
abilities for the two modes. The insets show details in the 
vicinity of the vacuum resonance. Note that the low to high 
cross section mode switching probability is enhanced around 
the vacuum resonance and reaches nearly unity. Because of 
this, scattering in the low mode can be dominated by mode 
switching, as it is in many of the results we describe in § 4. 

Absorption can also contribute significantly to the opac- 
ity over the frequencies and densities of interest, depending 
on the composition of the surface layers and magnetosphere. 
If the primary atom has a high atomic number, e.g., sili- 
con or iron, then even when kT — 20 keV the atoms will 
not be completely ionised (atomic binding energies are in- 
creased greatly by strong magnetic fields). If the primary 
atom is hydrogen or helium, virtually all atoms will be com- 
pletely ionised and only scattering and free-free absorption 
contribute to the opacity. Here we assume that only lighter 
elements are present, so that bound- free and bound-bound 
absorption may be ignored. 

Given this assumption, scattering dominates the opac- 
ity over most of the photosphere. Free-free absorption can, 
however, be important for the lower photon energies. The 
free-free absorption cross section in the low cross section 
mode is (see, e.g., Meszaros 1992) 



(Tabs w 4 x 10" 



ffi 



iV 26 T £ 



-1/2,. 



(64) 



where T 8 = T/10 8 K and lu w = hu/10 keV. When the 
photon energy hui < kT, induced scattering and absorption 
are important, and we include this by multiplying the total 
interaction cross section by the factor e n " / kT / ( e ntjJ / kT _ \y 

2.4 Summary of the effects of vacuum polarisation 

We have shown above that vacuum polarisation affects 
strongly the scattering cross section for the two modes. How- 
ever, vacuum polarisation does not introduce additional ab- 
sorption. As can be seen from equations (^), and (^), the po- 
larisation averaged cross section depends only on the imag- 
inary part of 72/ . From equation (J32J) , the vacuum contri- 
bution, represented by the third term in this equation, is 
real. Moreover, the imaginary part of ni is proportional to 
N e . Therefore, from equation (o), the polarisation averaged 
cross section depends on neither plasma density nor on vac- 
uum polarisation. The main effect of vacuum polarisation is 
redistribution of the cross section between the two modes. 
The effects described here will not be seen in calculations 
that use polarisation averaged cross sections. 

For example, the scattering cross section in the O- 
mode has no resonance at the cyclotron frequency when the 
vacu um effects are neglected (Ventura, Nagel, & Meszaros 
197E). Since radiative transfer takes place predominantly in 
the low cross section mode, neglecting vacuum polarisation 
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would inhibit formation of cyclotron lines in the polarisation 
averaged spectra. 

Vacuum polarisation alters the mode switching proba- 
bilities, which also influences radiative transfer. In the case 
considered here, the increase of P2-.1 to almost unity in the 
wings of the resonance also increases the number of scatter- 
ings that a photon undergoes before leaving. This increases 
the chances of absorption, and enhances the effects of Comp- 
tonisation. 



3 QUALITATIVE TREATMENT OF 
SCATTERING WITH THE SECOND 
VACUUM FREQUENCY 

From § 2 it is clear that the scattering cross section depends 
strongly on the parameter 

sin 2 6 lob 1 



2 cos 9 lo (1-2/) 



( •" ( '• , 



V2(h)) 



(65) 



When |6| <C 1 the cross section is ~ <tt in either polar- 
isation mode, whereas when |6| 3> 1 one mode has a <C <tt- 
Because lo C 2 ~ AT 1 / 2 , as radiation propagates from regions 
of high density to regions of low density the resonance fre- 
quency changes. In effect, the resonance acts as a sliding win- 
dow of high opacity. Qualitatively, we expect that the reso- 
nance will have an impact on the flux at a given frequency 
if a) the optical depth through the resonance is greater than 
unity, and b) the optical depth to the surface without the 
resonance is less than unity. These requirements give us, re- 
spectively, the minimum and maximum photon energies at 
which the vacuum resonance changes the spectrum. 

We can estimate the fractional width e of the resonance 
by arbitrarily setting the bound between \b\ <C 1 and |6| 2> 1 
at |6| = 1 and demanding that at lo — lo C 2(1 ± e), |b| = 1. 
This gives 



= 0.02/rV ^-(l-2/) , (66) 
sin 



lob sin 2 1 

which is identical to equation (|46|). If e <C 1 and the atmo- 
sphere is exponential with scale height I, then the thickness 
of the resonance is « Ael. Assuming that the resonance cross 
section is a « <7t sin 2 8, if the magnetic field is parallel to the 
surface normal then the optical depth through the resonance 
is r = N(4e£)a T sin 2 0/cos0: 

r = 6WrW? (l-2/) , B < B c 



30£iw? (l-2/) 



B > B c 



(67) 
(68) 



where l\ = (1/10 cm). This gives the lower bound to the 
feature. To estimate the photon energy at which the optical 
depth without the resonance is unity, we solve 



n <jt 



sin 2 (9 



cos8 I \lob 



(- 



1 , 



(69) 



which gives Tilu « 30^ 1/4 keV when B < B c and tito w 
20£~ 1/4 h 1/4 keV when B > B c . 

To summarise, when the magnetic field is parallel to the 
surface normal then when B <C B c this simple analytical 
treatment gives lower and upper bounds uii and u)h of 



hLo t « 6h~ 1/3 £~ 1/3 (l - 2/)~ 1/3 keV 
huj h w 3O^~ 1/4 (sin0) _1 keV , 
whereas for B 3> B c 
Tiloi w 3£ 1 " 1/3 (1 - 2/)~ 1/3 keV , 
fw h « 2O^ 1 ~ 1/ V /4 (sin0)~ 1 keV . 



(70) 
(71) 

(72) 
(73) 



The scalings with magnetic field, scale height, and / are 
confirmed by our numerical results: typically Tiloi ~ 5 keV 
and htOh ~ 40 keV, if the radiation is produced very deep in 
the atmosphere. 



4 NUMERICAL METHOD AND RESULTS 
4.1 Method 

Following the procedure outlined in § 2.3 to generate cross 
sections and adding the corrections mentioned there for free- 
free absorption and induced Compton scattering, we use an 
adaptive mesh to produce a table of physical depth as a 
function of optical depth for various values of the magnetic 
field and scale height. The physical depth depends on both 
the frequency and the direction of the photon, and in partic- 
ular for angles nearly perpendicular to the direction of the 
magnetic field the vacuum resonance feature is very sharp, 
and thus must be integrated carefully. We achieve sufficient 
accuracy with 130 logarithmically spaced zones in photon 
frequency from Tllo — 10 _2,6 m e c 2 to Tllo — rn e c 2 , 50 linearly 
spaced zones in angle, and optical depths from 10 _1 to 10 4 , 
logarithmically spaced. 

To produce a spectrum, we use a Monte Carlo code 
to pass 100,000 photons through the atmosphere, starting 
each photon at a single layer deep in the atmosphere. For 
each scattering we generate an optical depth traveled, which 
combined with the direction and with the table of physical 
depth versus optical depth gives the new physical depth of 
the photon. If the photon escapes from the atmosphere we 
store its energy and propagation direction. Otherwise, we 
use the differential cross section to select its new energy 
and angle and determine if the photon switches polarisation 
modes, and repeat the process. 

It is essential to include the effects of Comptonisation 
and mode switching in these simulations. At temperatures 
of ~10 keV the energy change of a photon in a single scatter- 
ing is much greater than the width of the vacuum resonance. 
Thus, in the >50% of scatterings in the resonance when the 
photon switches modes, the photon will be Comptonised out 
of the resonance and may require several thousand scatter- 
ings to return to the low cross section mode, during which 
its fractional energy change can be of order unity and in the 
course of which the photon is likely to be absorbed. We ac- 
count for the effects of mode switching by checking at each 
scatter if the photon switches modes and, if it does, simulat- 
ing an absorption by drawing the new energy of the photon 
from a thermal distribution. We do not change the physical 
depth of the photon in this case, because the mean free path 
in the parallel mode is small enough that the net distance 
traveled in the parallel mode is negligible. To speed up the 
program we also limit the total number of scatterings per 
photon to 1000; in the results reported here, between 97% 
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Figure 2. Model spectra, as seen at the stellar surface (solid lines), showing the effects of varying photon production depth. The input 
spectrum in each case is a 10 keV blackbody (dotted lines), and the magnetic field is perpendicular to the surface. The magnetic field 
for all four panels is B = 10 B c , the scale height is Z=10 cm, and the pair fraction is / = 0. In panels a, b, c, and d, respectively, the 
photons are produced at 10, 11, 12, and 13 scale heights, where the number density at the surface is 6 X 10 21 cm~ 3 (corresponding to 
10~ 2 g cm -3 for a fully ionised hydrogen plasma). Note that the low-energy bound of the feature remains roughly constant, whereas the 
high-energy bound increases with increasing production depth. 



and 100% of the photons escape in fewer than 1000 scatter- 
ings. 

4.2 Results 

In this section we show spectra for several representative 
combinations of physical parameters. In every case the in- 
put spectrum is a blackbody at 10 keV. Figure 2 shows a 



sequence where the production depth of the photons (in 
units of the scale height) changes but the magnetic field 
is fixed at B = 10 B c and the scale height is 10 cm. Note 
that the high-energy falloff moves to higher energies with 
greater depth, whereas the low-energy edge of the feature is 
roughly constant in energy, as predicted by our treatment in 
§ 3. Figure 3 shows spectra for different magnetic fields, with 
B = 0.3, 0.5, 1.0, and 10 B c . In these the vacuum resonance 
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Figure 3. Model spectra, as seen at the stellar surface (solid lines), showing the effects of different magnetic fields. As in Figure 2, the 
input (dotted lines) is a fO keV blackbody, the magnetic field is perpendicular to the surface, i=10 cm, and / = 0, and here the photons 
are assumed to be produced deep in the atmosphere. Panels a, b, c, and d, respectively, show spectra for magnetic fields of B = 0.3 B c , 
B = 0.5 B c , B = B c , and B = 10 B c . This figure shows that the absorption-like feature will only be detectable for B>B C . 



feature is only clearly visible for B>B C ; at lower fields the 
contrast in cross section in and out of the vacuum resonance 
is not enough to create a significant feature. In Figure 4 we 
isolate the effects of changing the scale height, as I ranges 
from 1 cm to 1000 cm. From these we see that the larger 
the scale height, the broader the feature, to the extent that 
the vacuum resonance may appear as only a falloff at low 
and high energies instead of an absorption-like feature. In 
Figure 5 we concentrate on the effects of different pair den- 
sities. In these graphs B — 10 B c , the scale height is 10 cm, 



and the production depth is 13 scale heights. As the frac- 
tion / of pairs (defined earlier) increases from to 0.49, the 
high-energy edge of the feature remains fixed while the low- 
energy edge moves up, narrowing the trough and making it 
shallower. Figure 6 shows our last energy spectrum, in which 
the magnetic field is in the surface instead of normal to it 
(as it was in Figures 2-5), and we see that the low-energy 
part of the spectrum is flatter than it is when the field is 
along the surface normal. 

In Figure 7 we show the angular distribution of emer- 
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Figure 4. Model spectra, as seen at the stellar surface (solid lines), showing the effects of different scale heights. Here B = 10 B c , / = 0, 
the input (dotted lines) is a 10 kcV blackbody produced deep in the atmosphere, and the magnetic field is perpendicular to the surface. 
Panels a, b, c, and d, respectively, show spectra for scale heights of / = 1 cm, I = 10 cm, I = 100 cm, and I = 1000 cm. Note that the 
vacuum feature moves to lower energies for larger scale heights. Thus, an instrument with a low-energy spectral cutoff of, e.g., ~ 5 keV 
(such as ICE) might see only a low-energy dropoff during a burst when the scale height is large, but could see the whole absorption-like 
feature during the afterglow when the scale height has dropped. 



gent radiation, integrated over all photon energies, for a 
magnetic field B — 10 B c along the surface normal, a scale 
height of 10 cm, and energy generated deep in the atmo- 
sphere. For comparison we also show the angular distribu- 
tion of emergent radiation for isotropic scattering. The ra- 
diation is strongly peaked towards the normal; part of this 
comes from geometric effects (see, e.g., Chandrasekhar 1960, 
§25), and part comes from the angular dependences of scat- 



tering in a strong magnetic field. For a magnetic field in the 
surface, the emission also depends on azimuthal angle, but 
qualitatively the angle-dependent scattering and geometric 
effects partially cancel and the emergent radiation is closer 
to isotropic than if only geometric effects operated. We also 
note that, regardless of the orientation of the magnetic field, 
the photon energy spectrum is essentially independent of 
viewing angle. 
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Figure 5. Model spectra, as seen at the stellar surface (solid lines), showing the effects of different pair densities, fn each panel 
B = 10 B c , I = 10 cm, and the input (dotted lines) is a 10 keV blackbody produced deep in the atmosphere, and the magnetic field is 
perpendicular to the surface. In panels a, b, c, and d, respectively, the pair fraction (defined as the ratio of the number of positrons to 
the total number of leptons) is / = 0, / = 0.4, / = 0.45, and / = 0.49. As derived in § 3, the pair fraction affects only the low energy 
range of the vacuum feature. 



5 DISCUSSION 

Many of the spectra in the preceding section show prominent 
absorption- like features due to the vacuum resonance. As 
discussed in §4.1, these features owe their existence in large 
part to the influence of Comptonisation and mode switch- 
ing, which imply that photons that interact in the resonance 
undergo frequency diffusion and spread their energy over the 
spectrum, creating a relative deficit at the original energy. 



In contrast, if Comptonisation is neglected the photons un- 
dergo primarily spatial diffusion, with minimal effects on the 
spectrum unless the atmosphere is of uniform density (see, 
e.g., Ventura et al. 1979). Thus, in treatments of this effect 
that do not include Comptonisation (such as the Bezchast- 
nov et al. 1996 radiative transfer approach) the vacuum res- 
onance will be smeared out by density gradients, whereas 
in our spectra the breadth and visibility of the absorption- 
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Figure 6. Model spectra, as seen at the stellar surface (solid 
line), for a magnetic field in the surface of the star. The low- 
energy spectrum is flatter than when the field is perpendicular to 
the surface. 
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Figure 7. Angular distribution of emergent radiation (solid line) 
for a magnetic field perpendicular to the surface, showing a strong 
peak along the surface normal (p, = 1). For comparison, we also 
show a simulation of the angular distribution for isotropic scat- 
tering (dotted line). 



like feature are partially because of the density gradients, as 
explained in §3. 

In the preceding section we examined the vacuum res- 
onance features produced for magnetic fields, scale heights, 
and pair fractions like those expected in SGRs. We now sum- 
marise the qualitative effects of each parameter, and com- 
pare our theoretical spectra with the spectra seen in SGR 
1806-20. 

The physical quantity of greatest interest is the mag- 
netic field, because above ~ B c many exotic processes play 
an important role. Our results show that if a vacuum reso- 
nance feature is seen, the field must be stronger than ~ B c . 
Since the location and strength of the feature depend only 
weakly on B, observation of the feature would only place a 
lower bound on the magnetic field; this would imply that 
the field of the SGR was at least as strong as the strongest 
magnetic field inferred for a pulsar, 2 x 10 13 G. 

If the density structure of the atmosphere can be rep- 
resented as exponential, then the scale height affects the 
qualitative nature of the observed spectrum. For large scale 
heights £>1000 cm, the feature is manifested as a dropoff 
at ~ 40 — 50 keV, and the rise at lower frequencies is be- 
low 5 keV. If instead ^<10 cm, both the low and high energy 
wings around the feature are in principle visible above 5 keV. 
These two regimes might be appropriate to a burst (for large 
scale heights) and to afterglow (for small scale heights). In 
particular, detailed time-resolved spectra should show this 
transition as the luminosity drops. 

The depth of energy production affects the lower bound 
to the vacuum resonance feature, and for a shallow enough 
production depth (say, nine scale heights when B = 10 B c ) 
the feature is not present. This implies that we may not ex- 
pect to see such features in accretion-powered pulsars, even 
if their fields are very strong; typical stopping depths are 
<50 g cm~ 2 (Miller, Salpeter, & Wasserman 1987), corre- 
sponding to less than nine scale heights. 

The pair fraction / does not affect the location of the 
upper bound to the vacuum resonance feature, but does 
change the location of the lower bound. If the plasma den- 
sity is dominated by pairs, the feature is narrowed and de- 
tectability requires that the source of photons be deeper 
than when the electron-baryon contribution dominates. If 
the magnetic field and scale height are known from other 
information, the centroid and width of the resonance fea- 
ture would then tell us the pair density and the depth where 
photons were generated. 

How do the spectra of the bursts from SGR 1806-20 fit 
into this picture? None of the spectra have obvious spectral 
features that clearly indicate a vacuum resonance. There is 
an indication (Fenimore et al. 1994) that there is a low- 
energy falloff below ~ 15 keV; this could be explained by a 
large scale height in our model, but many explanations not 
involving strong magnetic fields are equally viable. The X- 
ray Timing Explorer, with its large area (3000 cm 2 at 3 keV) 
and good spectral resolution (18% at 6 keV) will, if it catches 
a SGR burst, give us much-needed spectral information in 
the range where features such as the vacuum resonance are 
expected to appear; the High-Energy Transient Experiment 
may likewise provide useful information, particularly with 
its Wide Field X-ray Monitor. 

In conclusion, in magnetic fields B~>B C and plasma den- 
sities comparable to those expected in SGRs photon scatter- 
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ing is strongly affected by the vacuum resonance. Since scat- 
tering opacity is likely to dominate over absorption opacity, 
this is an important effect which must be included in cal- 
culations of radiative transfer through strongly magnetised 
atmospheres. 

Acknowledgements 

We thank Carlo Graziani, George Pavlov, and Victor 
Bezchastnov for valuable discussions. This work was sup- 
ported in part by NASA grant NAG 5-2868 and, through 
the Compton Fellowship Program, by NASA grant 5-2687. 



REFERENCES 

Adler, S. L. 1971, Ann. Phys., 67, 599 

Bezchastnov, V. G., Pavlov, G. G., Shibanov, Yu. A., & Zavlin, 
V. E. 1996, Proceedings of the Third Huntsville Workshop on 
Gamma-Ray Bursts, in press 

Chandrasekhar, S. 1960, Radiative Transfer (New York: Dover) 

Evans, W. D., et al. 1980, ApJ, 237, L7 

Fenimore, E. E., Laros, J. G., & Ulmcr, A. 1994, ApJ, 432, 742 
Gnedin, Yu. N., & Pavlov, G. G. 1974, JETP, 38, 903 
Gnedin, Yu. N., Pavlov, G. G., & Shibanov, Yu. A. 1978, JETP 
Lett., 27, 305 

Hurley, K., Sommer, M., Kouveliotou, C, Fishman, G., Mecgan, 

C, Cline, T., Boer, M., & Niel, M. 1994, ApJ, 431, L31 
Katz, J. I. 1982, ApJ, 260, 371 

. 1993, in Compton Gamma-Ray Observatory, 

ed. M. Friedlander, N. Gehrels, & D. J. Macomb (New York: 
AIP), 1090 

. 1994, ApJ, 422, 248 

Kouveliotou, C. ct al. 1993, Nature, 362, 728 

Kouveliotou, C. et al. 1994, Nature, 368, 125 

Kulkarni, S. R., & Frail, D. A. 1993, Nature, 365, 33 

Kulkarni, S. R., Frail, D. A., Kassim, N. E., Murakami, T., & 
Vasisht, G. 1994, Nature, 368, 129 

Lamb, D. Q. 1982, in Gamma-Ray Transients and Related Astro- 
physical Phenomena, cd. R. E. Ligenfelter, H. S. Hudson, & 

D. M. Worrall (New York: AIP), 249 

Mazets, E. P., & Golenetskii, S. V. 1981, Ap&SS, 75, 47 
Mazets, E. P., Golenetskii, S. V., Ilyinskii, V. N., Aptekar, R. L., 

& Guryan, Y. A. 1979, Nature, 282, 587 
Mcszaros, P. 1992, High-Energy Radiation from Magnetized Neu- 
tron Stars (Chicago: Univ. of Chicago Press) 

Miller, M. C, & Bulik, T. 1996, Proceedings of the Third 
Huntsville Workshop on Gamma-Ray Bursts, in press 

Miller, G. S., Salpeter, E. E., & Wasserman, I. 1987, ApJ, 314, 
215 

Pavlov, G. G., & Gnedin, Yu. N. 1984, Sov. Sci. Rev. E, 3, 197 
Pavlov, G. G., & Shibanov, Yu. A. 1979, ZhETF 76, 1457 
Pavlov, G. G., Shibanov, Yu. A., & Yakovlcv, D. 1980, Ap&SS, 
73, 33 

Thompson, C, & Duncan, R. C. 1995, MNRAS, 275, 255 

Tsai, W., & Erber, T. 1974, Phys. Rev., D10, 492 

Ventura, J., Nagel, W., & Meszaros, P. 1979, ApJ, 233, L125 



© 0000 RAS, MNRAS 000, 000-000 



